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Abstract. We determine the range of parameter space of an interacting quintessence (IQ) model 
that best fits the luminosity distance of type la supernovae data and the recent WMAP measurements 
of Cosmic Microwave Background temperature anisotropics. Models in which quintessence decays 
into dark matter provide a clean explanation for the coincidence problem. We focus on cosmological 
models of zero spatial curvature. We show that if the dark energy (DE) decays into cold dark matter 
(CDM) at a rate that brings the ratio of matter to dark energy constant at late times, the supernovae 
data are not sufficient to constrain the interaction parameter. On the contrary, WMAP data constrain 
it to be smaller than < 10^^ at the 3(7 level. Accurate measurements of the Hubble constant and 
the dark energy density, independent of the CMB data, would support/disprove this set of models. 

INTRODUCTION 

Recent observational data suggest the Universe has entered a period of accelerated 
expansion [1]. The ACDM model is the simplest model that best fits the supernovae 
luminosity distance and CMB anisotropy spectrum data [1, 2]. Though it looks a very 
serious candidate it raises two main theoretical problems. First, it assumes the energy 
density of the vacuum but according to quantum field theory it should be larger than 
observed by 120 orders of magnitude. Secondly, there is the coincidence problem, 
namely: why the energy densities of matter and dark energy -scaling so differently with 
expansion- happen to be of the same order precisely today? To solve the latter problem 
models featuring an interaction between quintessence and cold dark matter has been 
proposed [3, 4]. 



THE INTERACTING QUINTESSENCE MODEL 

To account for the coincidence problem Zimdahl et al. [3] postulated an interaction 
between cold dark matter and dark energy (a quintessence scalar field). 

p^ + 3H{l+w^)p^ = -5, p,d,n + ^Hpcdm = S. (1) 

This brings the ratio between their energy densities, r = pcdm/px, to a stable constant 
value at late times. The interaction term, 5, may depend on H, px and pcdm- A simple 
ansatz is 5 = 3Hc^{px + pcdm)^ where the constant is the interaction parameter. From 
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FIGURE 1. Evolution of the energy density ratio r = pcdm/Px from an unstable maximum towards a 
stable minimum (at late times) for different values of c^. As current value we take ro = 0.42. 



Eqs. (1) the aforesaid ratio is constrained to evolve between two stationary values, 

r± = — ± + where r_|_ denotes the unstable value and r_ the stable one 
-see Ref. [4] for details. As seen from Fig. 1 the stronger the interaction, the smaller 
Pcdm in the P^st. Obviously, such a dissimilar evolution of p^^^ as compared to the 
ACDM model, gives rise to a different expansion rate and a different evolution of the 
metric perturbations. The former modifies the luminosity distance of SNIa. The latter 
should be observed in the CMB anisotropy spectrum. Below we constrain the IQ model 
with the SNIa data of Riess et al. [1] and the CMB anisotropy data collected by the 
WMAP satellite [2]. 



OBSERVATIONAL CONSTRAINTS 

To constrain the IQ model with the SNIa data we assume the following priors: a spatially 
flat model (Qk = 0) and —1 < Wx < —0.6. As seen from Fig. 2 the lack of data for 
redshifts z > 1.8 does not constrain [5]. 

To constrain the model with the WMAP data we adopt a phenomenological point 
of view and impose that the interaction exists along the streamlines of the fluids (dark 
energy and CDM), i.e., 

U^y^T^^cdm) = -U^^nPcdm-^H{pcdm+Pcdm) = -S . (2) 

The net momentum transfer between CDM and dark energy alongside the accelerated 
expansion, limits the growth of the CDM perturbations. We implement the perturbations 
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FIGURE 2. Joint confidence intervals at 68%, 95 % and 99 .9% C.L. of the IQ model fitted to the "gold" 
sample of SNIa data of Riess et al. [1]. For convenience, the axis is represented using a logarithmic 
scale and it has been cut to < 10~^, though models with = have been included in the analysis. 



TABLE 1. Mean values and 1(7 C.L. of the parameters of the IQ and A 
CDM models. 



Parameter 


IQ model 


ACDM model * 


Interaction parameter, 


< 5 X 10-3+ 




CDM density/critical density, Q.c 


0.4±0.1 


0.29 ±0.07 


Baryon density/critical density, Q.i, 


0.063 ±0.012 


0.047 ±0.006 


Hubble constant, ho 


0.61 ±0.06 


0.72 ±0.05 


Spectral tilt, 


0.96 ±0.03 


0.99 ±0.04 


Amplitude, A 


0.9±0.1 


0.9±0.1 



* WMAP team results [2]. 

+ Upper limit at 1(T C.L. 

** This is a prior of the ACDM model. 



equations (2) together with V^Ty = —'^ fxTy i^cdm) iri the CMBfast code [6] and as- 
sume a flat model with {Q.]^ = 0), universes older than 12 Gyrs, BBN, and —1 <Wx< 0.6. 
By means of a Monte Carlo Markov chain we reconstruct the surface of the likelihood in 
the 7D-space parameter defined by the cosmological parameters c^, Q^x, Wx, Q.b,HQ, n^, A. 
Table 1 summarizes our results. Due to the interaction the CDM redshift is smaller than 
in the ACDM model, therefore a higher value of Q-cdm is needed to have enough gravita- 
tional potential at matter-radiation decoupling. Moreover, in order not to alter the height 
of the acoustic peaks the Q^i, value increases. Otherwise, the small mean value of Hq 
would suppress them. The likelihood shows clearly that the interaction parameter 
is compatible with the WMAP data (see Fig. 3). Nevertheless, non-interacting models 
are also allowed, unless the amplitude. A, is not a parameter (see Ref. [5]). Thus, is 
bounded from above. At la confidence level < 5 x 10^. The spectral tilt, Hs, and the 
amplitude of the perturbations at the horizon size. A, are in perfect agreement with those 
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FIGURE 3. Crosses represent the likelihood computed from 50,000 models, marginalized over all 
parameters but one. The solid line represents the best Gaussian fit. 



ofthe WMAPteam [2]. 



CONCLUSIONS 

We have presented an interacting cosmological model that solves the coincidence prob- 
lem. Despite the seemingly odd results of Tab. 1 it is in perfect agreement with SNIa data 
[1] and CMB data [2]. Moreover, the lack of power at the lower multipoles becomes a 
less serious problem. Indeed, because of the smaller redshift evolution of CDM energy 
density the gravitational potential evolves slower than in the ACDM model. Thus, a 
smaller Integrated Sachs-Wolfe effect produces less power in the CMB anisotropy spec- 
trum. A next generation of accurate measurements of the Hubble parameter and the dark 
energy density, independent of the CMB data, are needed to support/disprove this kind 
of models. 
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